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Recent neutron scattering experiments in the superconducting state of YBCO have been interpreted in terms 
of a magnetic collective mode whose dispersion relative to the commensurate wavevector has a curvature 
opposite in sign to a conventional magnon dispersion. The purpose of this article is to demonstrate that 
simple linear response calculations are in support of a collective mode interpretation, and to explain why 
the dispersion has the curvature it does. 

PACS numbers: 25.40.Fq, 71.f8.+y, 74.25.Jb, 74.72.-h 



Inelastic neutron scattering, which measures the mo- 
mentum and energy dependence of the dynamic spin 
susceptibility, has emerged as a powerful probe of high 
temperature cuprate superconductorad. One of the most 
striking features of the data is the appearance below the 
superconducting transition temperature of a commensu- 
rate resonance, in the case of the bilayer cuprates, YBCO 
and Bi2212B 13. In addition, an incommensurate mag- 
netic response has been seenJor energies below the res- 
onance energy in both casesQ~l3. Recently, it has been 
argued that these two structures are part of the same 
collective mode dispersionEl. This interpretation is un- 
usual in that the curvature of the purported collective 
mode dispersion relative to the commensurate wavevec- 
tor is opposite in sign to that of a conventional magnon 
dispersion. 

One of the motivations of this interpretation is the sim- 
ilarity of the measured dispersion to thjeorjetical results 
based on linear response theory (RPA)ll3il3. In these 
works, though, it was unclear whether the incommensu- 
rate response below resonance was actually due to a pole 
in the response function as would be expected for a col- 
lective mode interpretation. The purpose of this article 
is to investigate this issue in detail. 

The methodology is standard, and details can be found 
in Ref. An effective quasiparticle dispersion is used 
to calculate the bare bubble in the superconducting 
state, xo- The RPA susceptibility is then obtained as 
Xo/(l ~ "^Xo)- The most crucial input is the form of the 
dispersion. In Ref. |l^, two dispersions were analyzed, one 
based on angle resolved photoemission (ARPES) data in 
the normal state (tbl), another on ARPES data in the 
superconducting state where the Fermi surface was flat- 
tened around the node to enhance the effects of incom- 
mensurability (tb2). In addition, in this paper, a slight 
modification of this second dispersion was made to fur- 
ther enhance the incommensurability (tb3). All three dis- 
persions are listed in Table I for completeness sake. For 
the superconducting gap, the standard cos{kx) — cos{ky) 
form was assumed, with a maximum value. A,™,, taken 
from recent tunneling measurements in YBCOcJ. J was 
adjusted to obtain a resonance at 41 meV, the resulting 
J being similar in value to estimates based on neutron 



scattering in underdoped samplesEl. 

To begin, it is helpful to discuss the frequency depen- 
dence of the bare bubble, xo- In Figure 1, Rexo and 
Imxo are shown at a typical incommensurate wavevec- 
tor for the three dispersions of Table 1. In Imxo, there 
are two steps, a low energy one at threshold, and another 
at slightly higher energy. The two energies correspond to 
the spin gap energies Wg(±g) = minfc(i?fe -t- Ek±q), where 
k is confined to the first quadrant of the Brillouin zone 
and Ek are the quasiparticle energies in the supercon- 
ducting state. The steps are a consequence of the coher- 
ence factors associated with a d-wave order parameter. 
For the commensurate wavevector, Q = (tt, tt), these two 
energies are degenerate, so there is only one step. By 
the Kramers-Kronig relations, these steps lead to loga- 
rithmic divergences in Rexo which become simple peaks 
when damping is taken into account. (All plotted re- 
sults are obtained by replacing lu hy uj + iT in the energy 
denominators when calculating xo)- 

In Figure 2, Rexo and Imxo are shown as a function 
of q at an energy of 35 meV for the three dispersions of 
Table 1. In the first case (tbl), one finds the typical re- 
sult that Rexo is quadratic relative to the commensurate 
wavevector Q with a negative curvature (i.e., a maximum 
at Q). Since Rexo increases with uj (up to the spin gap 
energy ujg{Q)), one expects that the RPA pole condition, 
1 — JRexo = 0, will be satisfied first at Q, then at higher 
energies at wavevectors increasingly displaced from Q. 
This would then give a conventional magnon like disper- 
sion. The actual situation is more complicated, though. 
Note in Figure 2a that there is a second maximum in 
Rexo at an incommensurate wavevector, associated with 
i-Ug{q) as discussed in the context of Figure 1. As the en- 
ergy increases, this second maximum becomes the global 
maximum. 

In Figure 3a, the q vector where the RPA Imx is max- 
imal at a given lu (w(g)) is plotted for this dispersion. In 
agreement with experiment, one finds incommensurabil- 
ity below the resonance energy, but the RPA pole condi- 
tion is not satisfied below resonance. Instead the incom- 
mensurability simply tracks LUg(q). As a consequence, 
the intensity plummets below resonance, as can be seen 
in Figure 4a. Above resonance, the RPA pole condition 
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is satisfied, but as can be seen in Figure 3a, the disper- 
sion crosses into the continuum (i.e., beyond 0Jg{q)) just 
above resonance, so the mode becomes damped, as obvi- 
ous from the intensity plot of Figure 4a. 

The above behavior is not what is observed, in that 
the experimental intensity decays more slowly below res- 
onance. This is because the incommensurability effects 
associated with dispersion tbl are too weak. This can be 
contrasted with the dispersion analyzed by Brinkmann 
and LeellJ. In that case, the constant energy contours of 
Ek are quite flat for energies below l^max, leading to an 
enhancement in the incommensurate effects. This can 
be seen in Figures 2b and 2c, where Rexo is plotted 
for the other two dispersions considered here (tb2 and 
tb3). hirtiic first case (tb2), the dispersion used in ear- 
lier worko, one finds a rather flat behavior around Q due 
to the flat quasiparticle dispersion around k — (tt, 0). In 
addition, a global maximum now occurs at an incom- 
mensurate wavevector due to the flattening of the Fermi 
surface around the d-wave nodal (tt, tt) directions relative 
to dispersion tbl. 

The resulting dispersion, Lo{q)^ of Imx from tb2 is plot- 
ted in Figure 3b. This is inJaetter agreement with ex- 
perimental results in YBCOeIIj. In addition, it can be 
seen from this plot that the RPA pole condition is now 
satisfied below the resonance energy for a few meV until 
the spin gap energy, ujg{q)^ is encountered. For energies 
below this, the mode is strongly diminished in intensity 
as it traces out the edge of the continuum. This result is 
consistent with earlier workES. 

Although Figure 3b is in good agreement with experi- 
ment as far as w(g) is concerned, it is still deficient in 
that the intensity of Imx below resonance still drops 
off faster than exjperiment, as plotted in Figure 4b and 
also noted earlierliJ. To analyze this further, dispersion 
tb2 was modified by making the dispersion less flat in 
one direction around (tt, 0), and flattening the Fermi sur- 
face around the d-wave node even more. The resulting 
ReXo is shown in Figure 2c, and is siiiiilar to that ob- 
tained earlier by Brinkmann and Ledl3. One now finds 
a truly quadratic behavior around Q with positive cur- 
vature, a direct consequence of the stronger incommen- 
surate peaks. As a consequence, the RPA pole condition 
is satisfied for an even greater energy range below res- 
onance (Figure 3c), and moreover, the intensity below 
resonance falls off much more slowly (Figure 4c), as ob- 
served experimentally. 

The results in this paper concentrated on the (tt, qy) 
direction, since this is where the maximum in the in- 
commensurability was observed. Two dimensional q 
plots from -dispersions tb2 and tb3 are consistent with 
experimentD in showing a baseball diamond shaped in- 
commensurability pattern below resonance with global 
maxima along (tt, g^) and ((72,,7r), as reported earlierEJ. 
Recently, an anisotropy was observed in thiSpintensity 
pattern for partially detwinned YBCO samplestH. If this 
is a true magnetic anisotropy not related to the phonon 
background, then within the RPA context, it would have 



to be due to the influence of the chain bands wiping 
out two of the incommensurate spots, not surprising 
given the sensitivity of the incommensurability to the 
electronic structure. On the other hand, this is conjec- 
ture at this stage, as a proper calculation taking into 
account both plane and chain bands has yet to be per- 
formed. Of course, the most natural explanation of this 
anisotropy would be due to stripe formation, but it is 
unclear whether such a picture can reproduce the results 
demonstrated here (in particular, the energy dependence 
of the incommensurability). 

In conclusion, RPA calculations of the dynamic sus- 
ceptibility are in support of a collective mode interpreta- 
tion of the magnetic dispersion relation observed in bi- 
layer cuprate superconductors. The observed incommen- 
surability has strong implications for the effective quasi- 
particle dispersion, and it would be of interest to verify 
whether the dispersions analyzed here, in particular in 
regards to the flattening of the Fermi surface around the 
d-wave node, are consistent with high resolution angle 
resolved photoemission data, particularly in the case of 
YBCO. 

The author thanks Andrei Chubukov and Oleg Tch- 
ernyshyov for discussions. This work was supported by 
the U.S. Dept. of Energy, Basic Energy Sciences, under 
Contract No. W-31-109-ENG-38. 
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FIG. 1. -Rexo and Imxo for q = (l,0.8)7r versus ui for 
the three dispersions of Table I, with Amax~'29 meV (a and 
b) and 25 meV (c), r=0.5 meV, and T^ISK. Units are 
states/eV/CuO plane, and so should be multiplied by 2fi% 
to compare to neutron susceptibilities. 
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FIG. 3. q = {l,qy)Tr where Imx is maximal at a given 
uj (solid lines). The dashed lines are the spin gap energies, 
ujg{±q), defined in the text, and the open diamonds represent 
the RPA pole condition, 1 — JRexo{q) ~ 0. Same notation as 
Fig. 1, with J=262 meV (a), 111 meV (b), and 155 meV (c). 
Note that Imxo ~ for u> below the lowest dashed line. 
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FIG. 4. Value of Imx along the uj{q) dispersion relation 
from Fig. 3. 



TABLE I. Tight binding dispersions based on angle re- 
solved photoemission data. The first three columns list the 
coefficient, a, of each term (eV), that is e{k) =r^Cirii{k), 
with tbl and tb2 previously considered dispersiondij^ and tb3 
a modified version of tb2 as discussed in the text. The last 
column lists the basis functions (the lattice constant a is set 
to unity). 



tbl 



tb2 



tb3 



rii{k) 



0.1305 0.0879 

-0.5951 -0.5547 

0.1636 0.1327 

-0.0519 0.0132 



0.1197 
-0.5881 
0.1461 
0.0095 



^{coskx + cos ky) 
cos kx cos ky 
I (cos 2kx + cos 2ky) 
-0.1117 -0.1849 -0.1298 i(cos 2fca: cos fc;, -I- cos fc^^, cos 2fc;,) 
0.0510 0.0265 0.0069 cos 2fc^ cos 2fc„ 
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FIG. 2. Rexo and Imxo for uj—35 meV versus q = (1, qy)n. 
Same notation as Fig. 1. 
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